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ABSTRACT 

The Rossi X-ray Timing Explorer (RXTE) has made feasible for the first time the search for high-frequency (ft 
100 Hz) periodic features in black hole candidate (BHC) systems. Such a feature, with a 67 Hz frequency, recently 
has been discovered in the BHC GRS 1915+105 (Morgan, Remillard, & Greiner). This feature is weak (rms 
variability ~ 0.3% — 1.6%), stable in frequency (to within ~ 2 Hz) despite appreciable luminosity fluctuations, 
and narrow (quality factor Q ~ 20). Several of these properties are what one expects for a "diskoseismic" 
g— mode in an accretion disk about a 10.6 M© (nonrotating)— 36.3 M Q (maximally rotating) black hole (if we are 
observing the fundamental mode frequency). We explore this possibility by considering the expected luminosity 
modulation, as well as possible excitation and growth mechanisms — including turbulent excitation, damping, and 
"negative" radiation damping. We conclude that a diskoseismic interpretation of the observations is viable. 

Subject headings: black hole physics — X-rays: Stars 



1. INTRODUCTION 

In a series of previous papers (Nowak & Wagoner 1992; 
Nowak & Wagoner 1993; Perez et al. 1997; hereafter NW92, 
NW93, P97, respectively), we have discussed a class of modes 
in thin, Keplerian accretion disks that only exist in the pres- 
ence of strong-field gravity (not in Newtonian gravity), some of 
whose eigenfrequencies depend primarily upon the mass and 
angular momentum of the black hole. The modes are trapped 
by general-relativistic modification of the radial (k) and ver- 
tical (fljj epicyclic frequencies of free particle, circular orbit 
perturbations (P97, and references therein). 

These modes of oscillation are perturbations of the disk 
which are proportional to exp[i(at + m(j))]. With disk angular 
velocity £!(r), their corotating frequency is ui(r) = a + mQ. 
Three classes of modes (p— , g—, and c— modes) have been 
identified (P97, and references therein). However, here we shall 
be concerned solely with radial (m = 0) g— modes, which are 
the most relevant observationally 4 . These g— modes are trapped 
where uj 2 < k 2 , in the region where n achieves its maximum 
value (at r — 8 GM/c 2 for a nonrotating black hole) because, 
unlike in Newtonian gravity, «(r) rolls over at small r and van- 
ishes at the inner disk edge, thereby creating a resonant cavity. 
The modes with the fewest radial nodes have relatively large 
radial extents A?' w GM/c 2 . 

The frequencies (/ = —a/2ir) of the radial g-modes are 
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disk enter only through the small correction term involving the 
disk thickness 2h, and the radial (n) and vertical (j) mode num- 
bers, with 6 ~ 1 (as derived from the WKB solutions of P97). 
For thin disk models where h/r ~ 0.1 L/LEdd (in the mode 
trapping region), e n j is typically on the order of a percent for 
the lowest modes (n ~ j ~ 1). Therefore, the mode frequency 
is relatively independent of disk luminosity. (The mode width, 
however, scales as Ar oc y/cZ oc \/h, where c s is the sound 
speed in the disk.) The radial g— mode is thus a good candidate 
for explaining the recently discovered high-frequency features 
in the black hole candidate GRS 1 9 1 5+ 1 05 . 

The observed frequency of this feature is 67 Hz, with a full 
width half maximum of 3 — 4 Hz. (Morgan, Remillard, & 
Greiner 1996, 1997). This implies an effective Q = f/Af ~ 
20. Despite factors of ~ 2 luminosity variations in the source, 
the frequency remained constant to within 1 2 Hz. The root 
mean square (rms) variability of the 67 Hz feature varied from 
0.3% — 1.6% of the total observed X-ray luminosity, with the 
lower limit essentially being the Poisson noise limit. At its 
strongest, the rms variability ~ 1% below 6 keV, and the rms 
variability ~ 6% in the 12 — 25 keV band, indicating that the 
spectrum of the feature is harder than the integrated source 
spectrum. (The spectrum of the source is fit by a power-law 
exponentially cut-off at ~ 5 keV.) 

Below we discuss to what extent g— modes can explain these 
observations. In §2 we further review the properties of these 
modes, and make simple estimates of their luminosity modula- 
tion. In §3 we discuss application to GRS 1915+105. In §4 we 
discuss the role of turbulence and radiation as possible damp- 
ing, growth, or excitation mechanisms. In §5 we present our 
conclusions. 



(cf. P97). Here F(a) is a known function of the dimensionless 
black hole angular momentum parameter a = cJ/GM 2 , rang- 
ing from F(0) = 1 to F(0.998) = 3.443. The properties of the 
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4 These modes cover the largest area of the disk in a region where the disk temperature is maximum (away from the uncertain conditions at the disk edge). 
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2. LUMINOSITY MODULATION OF MODES 



For both the psuedo-Newtonian (NW92, NW93) and fully 
relativistic (P97) diskoseismology calculations, the equations 
are formulated in terms of the potential V = SP/p, where 
SP is the Eulerian variation of the pressure and p is the un- 
perturbed density. For modes with an azimuthal and time de- 
pendence oc exp[i(m</> + at)], the potential V(r, r/, <j), t) can be 
separated into a radial, V r (r), and a vertical, V v (r]), component, 
where 77 = z/h(r) is the dimensionless vertical coordinate. The 
(small) fluid displacement vector, £, can be related to the Eule- 
rian potential V(r, r/, 4>, t) via the equations 

dV 

(uj 2 - k 2 ) -1 exp[i(m(j) + at)]V v - ' 



6(r, <f>,z,t) 
€z(r,<l>,z,t) 



(u 2 - N 2 )- 1 exjp[i{m4> + at)} V r 



dr 
dVn 
dz 



(2) 



where N z (77) is the vertical buoyancy frequency, which we shall 
set to zero (cf. P97 for a discussion of the role of non-vanishing 
bouyancy frequency). We shall usually use A to denote a La- 
grangian variation of a given quantity, where we follow the def- 
initions of NW92, NW93, and P97. (For m = 0, g-modes 
have uj 2 — a 2 ^ n 2 .) 

The radial component of the fluid displacement approxi- 
mately satisfies the WKB relation 

" 2 c 2 s (r, 0)^ + - * 2 )(^ 2 - Tfii) W = , (3) 
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(We have also taken the unperturbed vertical barotropic struc- 
ture to be in hydrostatic equilibrium.) 

For most simple a-models, the energy generation rate per 
unit volume is approximately aP(r,z)Cl(r), The modes not 
only perturb the pressure (we ignore possible perturbations to 
a, and one can show that perturbations to SI are negligible), 
but they also perturb the locations of the disk boundaries. The 
variation of the luminosity is therefore 
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where T(r) is a slowly varying separation function (akin to 
a separation constant in the WKB limit), c s (r, 0) is the speed 
of sound at the disk midplane, W = (k 2 — w 2 )~ 1 dV r /dr, ( ^ 
and f6 = for a nonrotating black hole. We have taken 
the perturbations to be adiabatic, which means that the La- 
grangian perturbation of the pressure and density are given by 
AP/P = -yAp/p w -7V • £. Utilizing the radial WKB ap- 
proximation employed in equation (||), the (approximate) com- 
ponents of the divergence become 

~ [(Tf! 2 -Lo 2 )/h 2 tt 2 uo 2 } exp(iat) V r (r)V v {ii) 



where we have employed the above expression for the La- 
grangian pressure change, and employed the WKB approxima- 
tion in each harmonic term above. 

The luminosity fluctuation in general will have contributions 
from the fundamental frequency, as well as the first two har- 
monics. We have calculated the luminosity modulation for 
modes in a disk with radial sound speed profile equal to that 
of a radiation-pressure dominated Shakura-Sunyaev disk with 
L/L,Edd = 0.3. The exact radial profile does not greatly effect 
the estimates, as we can replace the quantity / P(r, 0)Cl(r)dz 
with F(r), the disk energy flux, which is known from energy 
conservation. As the mode width Ar oc V~h, the magnitude of h 
does affect the luminosity estimates. If we define the maximum 
magnitude (achieved at a radius, r m ) of the vertical displace- 
ment vector 6( r m: z i t) t0 be = Hh(r m ), then for the case of 
the mode with one radial node in V r and two vertical nodes in 
Vrj {dZ, z /dri even about the midplane) 5 we have 



0.6% H sm(at + 6) - 9.8% H 2 sin 2 (>/j + (9) 
+ 0.2% H 3 sin 3 (at + 6) . (7) 



(Here we have used the psuedo-Newtonian approximation of 
NW92, NW93 for our calculations.) This is the modulation of 
the bolometric luminosity. As the mode exists in the inner, hot- 
ter regions of the accretion disk, the modulation in restricted 
high energy bandpasses will be greater. However, we have not 
properly accounted for radiative transfer effects. Dispersion in 
photon diffusion times over the extent of the mode will decrease 
the modulated luminosity for large optical depths. 
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3. APPLICATION TO GRS 1915+105 

Although the above luminosity modulation is not large, Ti ~ 
0.6 — 1.6 produces, at the fundamental frequency, rms variabil- 
ity - 0.3 - 1.6% [ = (0.6%W + 0.2%ft 3 )/%/2 ], as is seen 
for the feature in GRS 1915+105. The first harmonic can pro- 
duce this observed rms variability for H ~ 0.3 — 0.7. (The 
rms variability is the coefficient of the sin 2 term divided by 
V8.) The modulation seen in GRS 1915+105 increases to 6% 
if the energy band is restricted to the 12 — 25 keV range. This 
qualitatively agrees with what we expect for our modes. If we 
compare the bolometric luminosity variation of the mode to the 
bolometric luminosity from the region r w 6 — 20 GM/c 2 , the 
relative fractional variation increases to ~ 6%. 

For the particular radial and vertical disk structure assumed 
here, the rms variability in the first harmonic is more signifi- 
cant than the rms variability in the fundamental frequency. As 

5 The first "even" mode about the midplane is not trapped in the region of the epicyclic frequency maximum. The first "odd" mode leads to a vanishing vertical 
integral in the above estimate of luminosity modulation. 



aflP(r, z) dz 



(5) 



where P'(r',z') = P(r,z) + AP(r,z), r' = r + £ r (r), 
z' = z + 6( z )' an d r /' r Oi z o are the disk boundaries (ro can 
be taken to go to 00 without loss of generality). Transforming 
variables, the change in luminosity can be written as 
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so little is known, both observationally and theoretically, about 
disk vertical profiles we cannot definitively state whether this 
feature is generic. Based upon energy arguments, however, we 
can say that generically disk modes can only produce rms vari- 
ability ~ (9(1%) if Ti ~ 1. We regard the question of which 
frequency, the fundamental or first harmonic, to associate with 
the observed 67 Hz feature as an open issue. 

Equation ([[]) indicates that a 67 Hz (fundamental) modulation 
can be produced in a disk around a black hole with a mass of 
10.6 Mq for a nonrotating hole, and 36.3 M@ for a maximally 
rotating hole. There are few "natural" frequencies to invoke in 
a black hole system. If one were to appeal to a process occuring 
at the marginally stable orbit (r = 6 GM / c 2 for a nonrotating 
black hole), not only would the luminosity from this region be 
lower than our mode estimates (due to the no-torque bound- 
ary condition used in disk models), but also the required hole 
mass would increase by a factor of at least 16/ (3v3) ~ 3 (if 
the 67 Hz feature represents the fundamental mode frequency). 
Our nonrotating hole, fundamental frequency, estimate of the 
required mass would yield an observed X-ray luminosity (Mor- 
gan et al. 1997) <; 30% of the Eddington luminosity, consistent 
with our assumptions above. 

4. MODE EXCITATION AND SELECTION EFFECTS 

4. 1 . Turbulent Damping and Excitation 

It is possible to use a parametrized stress tensor to esti- 
mate the effects of turbulent viscosity on the modes (NW92, 
NW93). The canonical energy of a radial mode is E c ~ 
<J 2 p{^ 2 + £ 2 )AV, where AV is the volume occupied by the 
mode. Isotropic turbulence produces a rate of change dE c /dt = 
-E c /t, with t ~ \ aa [h 2 /\ 2 r + h 2 /X 2 z ] I" 1 and X r , X z , re- 
spectively, being the radial and vertical mode wavelengths. The 
corresponding quality factor is given by 



against the turbulent damping (NW93). The modes are excited 
to an amplitude of |£ z | ~ a(h/ X r ) 3 ' 2 h, and ~ a^J X r /h h, 
for isotropic and anisotropic viscosity, respectively (NW93). If 
turbulence is playing the dominant role in damping and excit- 
ing the modes, then we have the following constraints. For 
isotropic turbulence, the Q value is only large for a <C 1; 
however, this implies a correspondingly small amplitude. For 
anisotropic viscosity, not only can we tolerate a larger a, but 
we also achieve a larger mode amplitude, although achieving 
the observationally required nonlinear mode amplitude is diffi- 
cult in either case. 



4.2. Negative Radiative Damping 

As a first approximation, we took the modes to be adiabatic. 
In reality, we expect there to be small entropy changes due to 
various effects, the most notable one being radiative losses. If 
we have a radiation pressure dominated atmosphere, as is likely 
in high-luminosity disks, one properly should use 
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where s is the specific entropy. We can estimate 6 the effect of 
this term for a radiation pressure dominated atmosphere, where 
one has 
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(c is the speed of light, n es is the electron scattering opacity, 
and we have invoked vertical hydrostatic equilibrium.) 

For perturbations of scalars, the Lagrangian perturbation 
operator A commutes with total time derivatives D/Dt (cf. 
Lynden-Bell & Ostriker 1967). Using this fact, and the fact 
that for our modes DAs/Dt = iuiAs, we can show that 



as X z ~ h/ j and A r ~ yhr/n, where j and n are of order of 
the number of vertical and radial nodes in any particular eigen- 
function. Thus, for a <C 1, we can have high mode Q. 

The other contribution to the fractional width A/// of the 
corresponding feature in a power spectrum comes from the 
number of modes that are significantly excited. If we assume 
that this will include those whose wavelengths are greater than 
the maximum eddy size L ~ a x l 2 h divided by the correspond- 
ing mach number ~ a 1 / 2 , equation (jl|) gives 
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The minimum effective frequency width will be composed of 
this span of modes, each broadened by 1/Q ~ a. 

The above estimates are for isotropic viscosity. If the tur- 
bulence does not efficiently couple to the vertical gradients of 
the modes, then the mode Q value is increased by a factor 
~ (jX r /h) 2 (NW93). Aside from damping modes, turbulence 
can also potentially excite modes. Velocity perturbations in the 
disk, Sv, are made up of a mode component, 6vm, an d a turbu- 
lent component, 8vt- Viscous damping arises from terms of the 
form SvmiSvtj, while mode excitation arises from terms of the 
form SvTiSvTj (NW93). It is possible to make simple estimates 
of the magnitude of the turbulent excitation, and balance this 
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where r es is the electron scattering optical depth of the atmo- 
sphere. The term proportional to d\nil 2 /dr in the above is 
negligible in the WKB approximation. The Lagrangian varia- 
tion of the entropy is therefore directly proportional to the La- 
grangian variation of the pressure. 

We do not know the true radial and vertical structure of a 
disk model that correctly describes the spectral observations of 
GRS 1915+105. However, for a "standard" Shakura-Sunyaev 
a— disk model, one can show, with uo ~ k, that 0JT es h/c ~ 
a" 1 , and hence is likely to be 3> 1. Taking this to be the 
case, we can combine equation ( |l0| ) and equation (|l2|) to yield 
AP/P = j'Ap/p, where 



7 ' „ 7 (i + i _J^) .7(1 + -') 



(13) 



We have subsumed our ignorance of the disk's vertical structure 
into the parameter a' which we expect to be of 0(a). 



Throughout this section our estimates shall be based upon vertically averaged quantities rather than fully z— dependent quantities. 
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To estimate the effects that the non-adiabatic terms have upon 
our modes, we can substitute the above into equation (^) with 
d 2 W/dr 2 w —k 2 W and then expand it about the adiabatic so- 
lution. Specifically, we then have a dispersion relation that can 
be written in terms of a function G(ui) = 0, to which we are 
adding a function H(lo) due to radiation damping. If ojq is the 
unperturbed mode frequency and 5uj is its perturbation, we then 
have to first order 



dG 
duj 



Slu + H(lu ) = . 



The mode Q value then becomes 



8u> 



' H(lo ) 



2T / Q 



(14) 



(15) 



where on the right side of the above we have used equation (^]), 
and have taken the term Y£l± to be dominant (i.e. T > 1). 
The minus sign indicates mode growth; that is, the modes are 
unstable to radiative losses. As typically k^ 1 > h, T > 1, and 
a < 1, it is fairly easy to obtain \Q\ ^> 1. 

4.3. Selection Effects 

As discussed in §2, the modes must approach nonlinearity in 
order to be observationally relevant. If isotropic turbulent ex- 
citation determines the mode amplitude, then there is a natural 
maximum vertical perturbation amplitude with H -C 1. If the 
turbulence is anisotropic, on the other hand, then it is possible 
to achieve H ~ 1, so long as a ~ \Jhj\ r £ 1. For this case 
we expect Q values ~ a~ 3 . Radiative excitation also naturally 
leads to Ti ~ 1. At some point, currently unknown, damping 
effects and nonlinearities must limit this radiative mode growth. 
As a zeroth-order approximation we take the linear growth rates 
calculated above to be relevant to the near nonlinear regime, 
and thus be equal to the saturated damping rate. The above 
radiative Q value then gives a rough estimate of the width of 
observationally relevant modes. 

The low rms variability of the modes means that only modes 
with high Q values will be detectable. Observational limits for 
Q ~ 20 were rms ~ 0.3%, and detectable rms oc Q -1 / 2 for 
narrow modes. Considering the dominant contributions as n 
and j increase, turbulence leads to Q oc j~ 2 (isotropic) and 
Q oc n~ 2 (anisotropic), whereas radiative damping leads to 
Q oc j 2 /n 2 The limiting detectable rms variabilty is therefore 
linear in n and j. Turbulent excitation only leads to high Q val- 
ues for modes with low j (isotropic) or low n (anisotropic). Ra- 
diative damping, on the other hand, only leads to high Q values 
for modes modes with low n and high j. We note, however, that 
the mode calculations require a smooth, well-behaved back- 
ground to perturb. If small scale turbulence truly is responsible 
for disk viscosity, our unperturbed hydrodynamic equations are 
only relevant for the largest scales, again favoring both low n 
and low j. For all cases considered above, a high Q value is 
most easily achieved for a -C 1. The 67 Hz mode seen in GRS 



1915+105, if a diskoseismic mode, therefore is likely one with 
low n and j, and a -C 1. 



5. CONCLUSIONS 

The main motivations for attributing the 67 Hz feature in 
GRS 1915+105 to a diskoseismic mode are that these modes: 
1) are related to a "natural" frequency in the disk (i.e. the max- 
imum epicyclic frequency); 2) their spectra are expected to be 
characteristic of the inner, hottest regions of the disk; 3) then- 
frequencies are relatively insensitive to changes in luminosity; 
and 4) they have low rms variability. This latter feature, al- 
though in agreement with the observations, is the strongest con- 
straint. These modes cannot be applied to systems that show 
^ 10% rms variability over a wide range of energy bands. We 
have identified two mechanisms, turbulent excitation and nega- 
tive radiation damping, that naturally lead to appreciable mode 
amplitudes with high Q value for a -C 1. 

Again, we stress that the only other "natural" disk frequency, 
the Keplerian rotation frequency at the last marginally stable 
circular orbit, would require a black hole mass a factor > 3 
larger. Furthermore, less flux is emitted from that region than 
from the g— mode region. 

We saw that there were two requirements for the g— modes 
to be observationally detectable. First, the g— modes had to be 
wide, or equivalently the disk had to be luminous, assuming that 
y/c^ oc \/L. Second, the mode amplitudes had to approach 
the nonlinear regime. Ideally, one should perform numerical 
simulations of the nonlinear development of the g— modes to 
determine mode widths, Q, and luminosity modulation more 
accurately. It is interesting to note that the MHD simulations of 
Stone et al. (1996) do show copious p— mode production asso- 
ciated with the turbulence (Gammie 1996). (MHD simulations 
have yet to be performed for rotation and epicyclic frequency 
profiles relevant to the g— mode trapping region.) 

These g— mode oscillations should also exist in accretion 
disks around compact (soft equation of state), weakly magne- 
tized (B < 10 8 gauss) neutron stars. Under these conditions, 
the inner radius of the disk will be less than those radii where 
the g— mode exists. Again, large amplitude (rms 3> 1%) fea- 
tures cannot be explained with these modes. 

Even if the 67 Hz feature seen in GRS 1915+105 does not 
turn out to be a diskoseismic mode, it points out two important 
lessons. First, BHC systems can produce relatively stable, high- 
frequency features. Second, the Rossi X-ray Timing Explorer is 
capable of detecting and characterizing these features despite 
their weak variability. The search for diskoseismic modes in 
this system and other BHC has therefore become a viable and 
worthwhile pursuit. 
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